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A new nonlocal form of the off-shell kinetic equation is derived. While being equiv- 
alent to the Kadanoff-Baym and Botermans-Malfliet formulations in the range of for- 
mal applicability, it has certain advantages beyond this range. It possesses more ac- 
curate conservation laws for Noether quantities than those in the Botermans-Malfliet 
formulation. At the same time the nonlocal form, similarly to the Botermans- 
Malfliet one, allows application of the test-particle method for its numerical solu- 
tion, which makes it practical for simulations of heavy-ion collisions. The physical 
meaning of the time-space nonlocality is clarified. 
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I. INTRODUCTION 



The work of S.T. Belyaev and G.I. Budker, who demonstrated the Lorentz invariance of 

;he relativistic distribution function and derived relativistic Fokker-PIanck kinetic equation 
1 

stands in the line of achievements of the kinetic theory. This work entered many text- 
books and found numerous applications in various fields. Presently the relativistic transport 
concepts are a conventional tool to analyze the dynamics of dense and highly excited matter 
produced in relativistic heavy-ion collisions. 

A great progress was also achieved in microscopic foundation of the kinetic theory. The 
appropriate frame for description of non-equilibrium processes within the real-time formalism 
of quantum-field theory was developed by Schwinger, Kadanoff, Baym and Keldysh 
^. The formalism allows extensions of the quantum kinetic picture beyond conventional 
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approximations (like the quasiparticle one). This is caused by the quest for dynamical 
treatment of broad resonances as well as stable particles which acquire a considerable mass 
width because of collisional broadening. The above mentioned applications request for 
development of approximate self-consistent schemes possessing conservation laws being at 

nndnnnn 

least approximately satisfied [5|, l6|, 17|, l8|, l9|, llO, lUl. Based on these schemes, numerical 
transport methods for treatment of the off-shell dynamics have been developed jl^ . [isl . [l^ . 

Two slightly different forms of the Kadanoff-Baym equations expanded up to first-order 
space-time gradients are now used: the proper Kadanoff-Baym (KB) form, as it follows right 
after the gradient expansion of exact KB equations [sl, and the Botermans-Malfliet (BM) 
one, as it follows after a modification of a Poisson-bracket term in the KB equation |l5 |. 
Both the KB and BM forms coincide in the first-order gradient approximation but differ 
in higher orders. Both forms have their advantages and disadvantages llQl. The KB form 
possesses exact conservation laws for the Noether current and the energy-momentum [9|, 110[ , 
however, it does not allow the efficient test-particle method to be applied to its numeric 

n n n 

solution. The BM form is very suitable for the test-particle method ^2|, 113|, IMj but only 
approximately conserves the Noether current and the energy-momentum. 

In this paper we would like to put forward a new, space-time nonlocal form of the quantum 
kinetic equation, which combines favorable features of the above mentioned KB and BM 
forms. In sect. [TTlwe start with brief review of properties of the KB and BM forms of the 
off-shell kinetic schemes. Technical details are deferred to the Appendix El In sect. Illll the 
new nonlocal form of the quantum kinetic equation is derived and physical meaning of the 
nonlocality is clarified. Possible applications in numerical transport schemes are discussed. 



II. KADANOFF BAYM AND BOTERMANS MALFLIET KINETICS 

In this section we summarize the formulation of the off-shell kinetic equations in the two 
different forms: in the KB form and in the BM form. We assume the reader is familiar 
with the real-time formulation of non-equilibrium many-body theory and use of the contour 
matrix notation, detailed in Appendix |X1 

Starting point of all considerations is the set of Kadanoff-Baym equations which express 
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the space-time changes of the Wigner transformed^ correlation function iG ^{X,p) in terms 
of the real-time contour convolution of the self-energy E with the Green function G. We 
give the kinetic equation in compact notation (cf. Eq. (1A4I) ) 

v^a^iG-+(X,p) = [E®G-G®S]-+ with v^ = ^G^\p), (1) 

where Gq ^(p) is the Fourier transform of the inverse free Green function 

— for relativistic bosons 

Go\p) = { (2) 
I Po — p^/(2m) for non-rel. fermions or bosons. 

For a complete definition, Eq. ([1]) has to be supplemented with further equations, e.g., for 
the retarded Green function together with the retarded relations ( 1A6I) . If a system under 
consideration is only slightly spatially inhomogeneous and slowly evolving in time, a good 
approximation is provided by an expansion up to first order in space-time gradients. Then 
the main problem to arrive at a proper kinetic equation consists in accurately disentangling 
a rather complicated r.h.s. of Eq. ([1]). 



A. <l>-derivable approximations 

In actual calculations one often uses approximations or truncation schemes to the exact 
non-equilibrium theory, where conservation laws (such as charge and energy-momentum 
conservations) and thermodynamic consistency of the transport theory are not evident. It 



was shown 
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that there exists a class of self-consistent approximations, called $ 
derivable approximations, which are conserving at the expectation value level, i.e. they 
provide true Noether currents and a conserved energy-momentum tensor, and at the same 
time thermodynamically consistent. In these schemes the self-energies are self-consistently 
generated from a functional ^[G] through the following variational procedure [3] 

6i^G] J 2 for real fields 
iS,fe(X,p) ^ I ^ for complex fields ' ^' ^ ^ ^ + ' ^ 

The functional specifies a truncation scheme. It consists of a set of properly chosen 
closed two-particle irreducible diagrams, where lines denote the self-consistent propagators 



^ The space-time variable is X = = (t, x), and the Fourier transformed variable is p = 
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G, while vertices are bare. The functional variation with respect to G diagrammatically 
implies an opening of a propagator line of 

B. Physical notation 

It is useful to eliminate the imaginary factors inherent in the standard Green function 
formulation and introduce quantities which are real and possess clear physical meaning. 
Thus instead of Green functions G^^{X,p) and self-energies with i,j G { — h} 

in the Wigner representation we use the kinetic notation of Refs. [3, [l^. We define the 
generalized distribution functions F and F in the 8-dimensional phase space, 

F{X,p) = (T)iG-+(X,p), F(X,p) = iG+-(X,p). (4) 

Here and below the upper sign corresponds to fermions, while the lower sign, to bosons, 

A{X, p) = -2Im G^{X, p) = F±F (5) 

is the spectral function, and G^ is the retarded propagator. The spectral function satisfies 
the sum rule 

( '^^A{X,p) = 1 for nonrelativ. particles, (6) 
Jo 27r 

PqA{X,p) = 1 for relativ. bosons, 

, 27r 

which follows from the canonical equal-time (anti) commutation relations for (fermionic) 
bosonic field operators. The gain and loss rates of the collision term are defined as 

r-(x,p) = Tis-+(x,p), r-*(x,p) = is+-(x,p) (7) 

with the damping width 

T{X,p) = -2ImS^(X,p) = r°"*(X,p) ± r'^(X,p), (8) 

where E'^ is the retarded self-energy. 

In terms of above kinetic notation, the gradient-expanded Kadanoff-Baym equations 
are reduced to equations for real quantities: for the real and imaginary parts of the re- 
tarded Green function G^, and for the phase-space occupation F. Note that the number of 
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equations for G^^ Green functions is four, which twice exceeds the number of unknown func- 
tions {F and G^). Before the gradient expansion all equations were completely consistent. 
However, after the gradient expansion their interrelation is no longer obvious. Necessary 
interrelations have been derived in Refs. 

The equations for the retarded propagator in the first-order gradient approximation can 
be immediately solved with the result jJQ 

r(x,p) 



G 



R 



1 



M{X,p)+iT{X,p)/2 



AiX,p) 
ReG^{X,p) 



M2(X,p) + r2(X,p)/4 

M{X,p) 
M2(X,p) + r2(X,p)/4 



(9) 



with the "mass" function 



M(X,p) = Go^(p)-ReS^(X,p). 



(10) 



Note that the algebraic solution (Q is valid within the first-order gradient approximation, 
i.e. only the terms 0{dx) are omitted. 

C. Kadanoff Baym form of kinetic equation 



In terms of above notation, the KB kinetic equation for F in the first-order gradient 
approximation takes the form, see ^2] for details. 



VF{X,p) - {r'',ReG^} = C{X,p). 



(11) 



We refer this as the kinetic equation in the KB-choice^. Here {...,...} denotes the four- 
dimensional Poisson bracket 



V 



The differential drift operator is defined as 



dpi" J dXi" dp^' 

note that Dip = {M,iIj} for arbitrary function ip. The collision term is defined as 



(12) 



c{x,p) = r-(x,p)F(x,p)-r-*(x,p)F(x,p). 



(13) 



^ If the system consists of different particle species, there is a set of coupled kinetic equations corresponding 
to each species. 



6 



In terms of a functional $ the explicit form of the collision term is 

^<-^-'"^«|^^<-^-'"-^^'^-^'' 

cf. Eq. 

Within the consistent gradient expansion C oc dx- If the diagrams for the self-energy 
contain internal vertices, it gives rise to non-local effects in the collision term 

C = Cioc + Cmem- (15) 

Here Cioc is a local part of C, where all these nonlocalities are disregarded, and Cmem oc 



dx is the nonlocal part expanded up to first gradient terms. Following j?] we call these 
effects the memory effects. Explicit form of the Cmem depends on the specific system under 
consideration. 

For the sake of clarity, below we do few simplifications. We do not consider the memory 
effects, implying that C = C^°'^. Besides, we confine ourselves to the case void of derivative 
coupling. Also we do not explicitly introduce the particle-specie label to avoid overcompli- 
cation of equations. 

The local part of the collision term is charge and energy-momentum conserving by itself 

Tr / I " I Ci°^ = 0. (16) 



(2vr) 



Here e denotes a charge (e.g., the baryon number), while Tr implies the sum over all possible 
internal degrees of freedom (like spin, isospin, etc.) and over possible particle species. 
The true Noether current of the charge e, 

/(X) = eTr I ^^v^F{X,p), (17) 

is exactly conserved 

d,nX) = 0, (18) 

which follows right from the operator expression for this quantity, cf. Ref. Q]. 

On the other hand, the integration of the KB equation ffTTl) permits us to derive another, 
"effective KB current" 



JKB-eff(^) = eTr 



dp 



Vf^ 1 F - Re G 



dpf^ J dpfj, 



(19) 



7 



which is exactly conserved. The equivalence of [T7\ ) and [T^) and thus the exact con- 
servation of the Noether current f|T7|) follows from the corresponding invariance of the $ 
functional, cf. Eq. (6.9) in Gj, which results in the consistency relation 



eTr / <! -d, 



{271 



+ C\=0. (20) 



Similar situation takes place for the energy-momentum conservation. Here we also rely on 
^-derivable approximations which provide the corresponding exact conservation. 

The conserving feature is especially important for devising numerical simulation codes 
based on this kinetic equation. Indeed, if a test-particle method is used, one should be sure 
that the number of test particles is exactly conserved rather than approximately. In the 
test-particle method the distribution function is represented by an ensemble of test particles 
as follows 

F(X, p)r^Yl ^^'^ - W) ^^'^ (P - ' (21) 

i 

where the i-sum runs over test particles. Then the VF term in Eq. (fTT!) just corresponds to 
the classical motion of these test particles subjected to forces inferred from ReS^, while the 
collision term C gives stochastic change of test-particle's momenta, when their trajectories 
"cross" . For a direct application of this method, however, there is a particular problem with 
the KB kinetic equation. The additional term, i.e. the Poisson-bracket term |r™,ReG'^}, 
spoils this simplistic picture, since derivatives acting on the distribution function F appear 
here only indirectly and thus cannot be included in the coUisionless propagation of test 
particles. This problem, of course, does not prevent a direct solution of the KB kinetic 
equation. For instance, one can apply well developed lattice methods, which are, however, 
much more complicated and time-consuming as compared to the test particle approach. 



D. Botermans Malfliet form of kinetic equation 

As can be seen from Eqs. (|4]), ([8]) and f|T3l) . the gain rate differs from FT /A only by 
corrections of the first order in the gradients 

= rF/A + C/A = rF/A + 0{dx), (22) 

since C ~ 0{dx)- This fact permits us to neglect the correction 0{dx), as in the kinetic 
equation it leads to terms of already second-order in the gradients. Upon substitution 
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-pm _ YF/A^ proposed by Botermans and Malfliet [15|], one arrives at the following form of 
the kinetic equation, see [?] for details, 



PF(X,p)-<!r-,ReG^ 



(23) 



which is still equivalent to the KB form within the first-order gradient approximation (all 
terms oc 0{d\) are now omitted). We call the so obtained Eq. fl23l) the kinetic equation in 
BM-choice. All KB-choice properties of Eq. (1231) within a ^-derivable approximation also 
transcribe to BM-choice through the substitution = TF/A in the consistency relation 
Eq. ([20]). 

The BM equation exactly conserves the following "effective BM current" 



JBM-efr(^) = eTr 
= eTr 



dV r 

(2^2 



v^'F{X,p) + ReJ: 



R 



OF 

dPn 



-ReG 



Rdj TF/A) 

dPn 



(24) 



Bf'F, 



where 



A 



9ReS^\ MOT 



dp^ 



r dp, 



(25) 



is the flow spectral function introduced in 7|]. The 0-component of the flow spectral function 
B^ obeys the same sum-rule as the spectral function A, cf. jl^. The effective current ([21]) 
differs from the true Noether current of Eq. f[T7]) and the effective KB current f[T9l) in 
terms of the order of 0{dx), provided a ^-derivable approximation is used for self-energies. 
Thus with the BM choice the conservation laws of the Noether current ([T7)l and the energy- 
momentum tensor are only approximately fulfilled. ^ ^ 

The effective BM-current ( l24l) was used by S. Leupold [14] as a basis for the construction 
of a test-particle ansatz for numerical solution of the nonrelativistic BM kinetic equation. To 
automatically fulfill the effective current conservation, the test-particle ansatz is introduced 
for the combination 



ItBoF{x,p) ~ ^^'^ - ^^'^ - 



(26) 



rather than for the distribution function itself. Note that the energy p^{t) of the test particle 
is an independent coordinate, not restricted by a mass-shell condition. W. Gassing and S. 



Juchem 



13| used this test-particle ansatz in the relativistic case. 
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The BM kinetic equation (123|) together with ansatz (l26ll for the distribution function 
result in the following set of equations for evolution of parameters of the test particles 
between collisions 

*■ = .o-a..ReEH-(M/r)9..r + ^-^^ + <^'/r'^-r) • 

P- = .,„-a,.ReEH-(M/r)a,.r (^'^^^^ + ('''/r)v..r) . (28) 

^' = .o-a..ReEH-(M/r)9..r + (*'/n9.r) . (29) 



These equations o ' 
a test particle fl4 | 



motion, in particular, give the time evolution of the mass term M, of 



dMi _ Mi dTi 

the origin of which can be traced back to the additional term jPF/A, ReG'^} in the BM 
equation (l23l) . Here Mj(t) = M[t, Xj(t); pj(t)] measures an "off-shellness" of the test 
particle, and Ti(t) = T[t, Xj(t); Ei(t), Pi(t)]. Equation of motion ( |30l) yields Mi = ajFj, where 
ai do not depend on time, and implies that once the width drops in time the particles are 
driven towards the on-shell mass, i.e. to M = 0. This clarifies the meaning of the additional 
term jPF/A, ReG^} in the off-shell BM transport: it provides the time evolution of the 
off-shellness. 

The problem of not exactly conserving Noether charges within the BM kinetics can be 
mitigated in certain cases. For instance, if we start from initial conditions defined for 
particles in vacuum (where S = and hence both charges coincide) and end the evolution 
also in a very dilute state, the Noether charge turns out to be conserved in the end of the 
evolution. However, this is not the case with heavy-ion collisions, where we start with two 
cold nuclei. 

III. NONLOCAL FORM OF OFF-SHELL KINETIC EQUATION 

A. Nonlocal kinetic equation and conservation laws 

Let us rewrite the KB kinetic equation ( fTTj) in the following non-local (NL) form 
VF{X,p) - |r^,ReG'^| = C^^ 
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Eq. (pTl) is the key equation of our work. The coUision term C^^ is expressed here in terms 
of shifted variables. Note that the specific 8-phase-space memory in C^^ should not be 
confused with the memory effects resulting from internal structure of self-energy diagrams 
with more than two vertices. In order to distinguish between these two effects we call the 
non- locality in Eq. ( l3TI) . as delays (positive or negative) in 8-phase-space rather than the 
memory. 

To verify that the NL kinetic equation Eq. ( l3Ti) is equivalent to the KB one, including 
second-order gradient terms, let us expand the collision term C^^ in small quantities 

SX>' ^ iSrM) - j^^, (32) 



Then 



= (l + |i,ReG^}) (c + i{C,ReG«}) 



= C + {C/A,ReG^} + 0{d^x)- (33) 

Here we have taken into account that C oc 0{dx)- The resulting Poisson bracket 
|C/y4, Re G^} combines with the Poisson bracket on the l.h.s., cf. Eq. ( |22l) . resulting 
in |r™,ReG'^}, as it stands in the KB kinetic equation flTT]) . Thus, indeed, Eq. flMl) coin- 
cides with the KB equation ffTTl) within the approximation including second-order terms in 
space-time gradients. The roughening of this approximation by neglecting the second-order 
terms (i.e., if one puts C^^ = C), returns us to the BM equation. 
The consistency condition related to Eq. fl3ip can be rewritten as 



- 8p, J + 0(4)1^0. (34) 

cf. Eq. (!20l) . Thus, the precision of conservations of the Noether baryon number and the 
energy-momentum in the NL equation is one order in space-time gradients higher than 
that in the BM equation. Indeed, the effective BM current differs from the Noether one 
in terms of the first-order gradients whereas the effective current that follows from the NL 
kinetic equation differs only in second-order gradients. The latter takes place because the 
NL equation is equivalent to the KB equation including O 



^d\) terms. Let us recollect that 



Noether currents are exactly conserved for the KB choice 



a: 
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Within the range of formal apphcabihty, the NL equation is equivalent to the KB and 
BM forms discussed above. However, in practical calculations one frequently needs to apply 
the kinetic approach beyond the scope of its applicability. For instance, such a problem 
arises when one describes the initial stage of heavy ion collisions. We do not know how good 
the quantum kinetic equations are beyond the region of their formal applicability, but we 
certainly wish to respect conservation laws. From this point of view the KB form is certainly 
preferable. As the conservations are exact, we can still use the gradient approximation, 
relying on a minor role of this rather short initial stage of heavy ion collisions in the total 
evolution of a system. Exact conservation laws allow us to keep control of numerical codes. 
However, as we have mentioned, the efficient test-particle method is not applicable for 
solution of the KB equation. At the same time, for the BM kinetic equation the test- 
particle method is already available l^,[l^, for the price that it deals with the approximately 
conserved Noether current. 

When the collision term C is large, the neglect of the C /A term in the r.h.s. of the relation 
( I22l) is a bad approximation. Please, recollect that this approximation is in the basis of the 
BM approximation. Therefore, the nonlocal form of the kinetic equation which we present 
here is certainly preferable. Conservations of the Noether current and energy-momentum 
are more accurate here than in the BM form, since they are closer to exact conservations of 
the KB form.^ At the same time the test-particle method is applicable to the NL equation. 
In the NL case the set of equations for evolution of parameters of the test particles between 
collisions is the same as in the BM case, see Eqs. fl27l ) -fl29 l) . The only difference with the 
BM case is that collisions of test particles occur with certain time (and space) delay (or 
advance) as compared with the instant of their closest approach to each other. Because 
of this circumstance the effective BM current fl2^ is no longer locally conserved. Indeed, 
integrating C^^ over 4-momentum and taking into account relation (l33l) . we arrive at 



S„(-Ti/7^?^)+OK). (35) 



(2jr)M ap,. 

Precisely this contribution makes the KB effective current f[T^ from the corresponding BM 



^ Deriving (|31[) . we only added higher (than second order) space-time derivatives to equation in KB choice 
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current (12^ . This nonconservation does not prevent us from application of the test-particle 
method. Since collisions are nonlocal in the time-space, we need only global conservation 
of the BM effective current, which in fact takes place 

IV / <i'X |^,C- ^ / i'xa, (-T. / Ifj^^) . Om ^ Om (36) 

Now this conservation is approximate, 0{dx) terms are dropped. Nevertheless, this precision 
is better than that resulting from the formal accuracy of the kinetic equations, where already 
0{dx) terms are dropped in the approximately conserving quantities. 

Actually we could modify the NL equation in such a way that it exactly fulfills global 
conservation of the effective charge, which may be important for the test-particle method. 
Indeed, we recognize that to the first-order gradients the pre-factor J = (l + |^,ReG''^}) 
in Eq. (!3T|) is the Jacobian J of the transformation from d^pd^X to the d'^pdX integration 
expressed in terms of shifted variables pf^ = -\- 5p^ and = — 6X^^. By adding 
terms SJ of higher (than 0{dx)) orders to J, we do not violate the nature of the NL 
approximation. Then 6 J can be tuned in such a way that J +6 J becomes the exact Jacobian 
of the transformation and hence the BM effective current turns out to be exactly globally 
conserved: 

However the explicit expression for 6 J oc 0{dx) is very cumbersome. In order to keep the 
formalism transparent, we avoid this extra complication. 



B. Time delays and causality 

The NL kinetic equation flMl) contains the delayed collision term and the BM particle 
drift term including the drag and the back flow contributions. The extra Poisson-bracket 
term |C/A,ReG'^} which in the KB equation had a poorly defined physical meaning is 
now hidden in the 8-phase-space delays of the collision term thus acquiring a new physical 
interpretation. The space-time part of the retardation/advance (132|) is expressed as 

SX^ = Ib^-Uv^-?^], (38) 



2 r V 9p, 

where we used Eq. (|9]) and the definition of the flow spectral function (l25l) . There are 
two contributions to the 6X^. The first term is proportional to the same function 
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that determines the trajectories of the test particles between colhsions, see Eqs. (1271) - (l29l) . 
The second term (f ^ — (9Re S^/Sp^) /V relates to collisions: l/F ~ \/{u(yp) in the gas 
approximation, where u is the particle velocity, p is the density of the medium and a is the 
cross-section. 

Note that the time delay in the collision term has no definite sign (see discussion in 
the next subsection). Therefore, one could naively suppose that NL kinetic equation vio- 
lates causality. Here the following comment is in order. The initial Dyson equations for 
non-equilibrium Green functions are causal. The apparent problem appears already at the 
moment when one chooses t = {ti + t2)/2 and At = t2 — ti as independent variables (where 
t2 > ti due to causality) to perform then the Wigner transformation. Although approx- 
imations have not yet been done, in new variables the Green function G^{t, At) contains 
an information from the future, see Fig. 19 of Ref. [l7|. The problem would completely 
disappear, if one had used ti and At = t2 — ti as independent variables. Thus, working in 
t = {ti + 12)/2 and At = t2 — ti variables we incorporate some information from the future. 
Otherwise our results would be incompatible with those obtained from the original Dyson 
equations. However, this fact does not contradict causality, since the Dyson equations are 
causal and approximations have not yet been done. The necessity to include the informa- 
tion from the future is clearly seen, e.g., from analysis of the ^-diagrams containing more 
than two vertices. The vertices depend on different time variables. These dependencies 
reveal themselves in the above mentioned memory effects in the gradient expanded kinetic 
equations. 

The apparent acausality disappears after the first-order gradient expansion, and it does 
not manifest in the KB and BM forms of the kinetic equations since all variables there 
depend locally on the single time t. However via corresponding Poisson bracketed terms these 
equations contain an information on the mentioned non-local effects, which were originally 
present. Constructing ( !3TI) we re-grouped some of these terms hiding them in the collision 
term at the price that the apparent problem again appeared. However, the NL kinetic 
equation can be read from the right to the left in such a way that the collision term in 
future is determined by the drift term at present. This way the NL equation is again 
completely causal. 

As we will see in the next section the new NL form of the kinetic equation is convenient to 
discuss the so called negative time delays which have been experimentally observed in semi- 
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conductors, see [3] and refs. therein. The quantum kinetic delay effects reffect fundamental 
physical limitations which follow from basic uncertainty principles for the time-energy and 
the coordinate-momentum. The finite duration (positive or negative) implies that the final 
state of the system is not entirely determined until the process is completed. The wave-like 
nature of excitations indicates that the interactions between waves are interference phenom- 
ena which have finite duration on a scale of at least one period of the interacting waves. 
Collisions cannot be considered to be completed at times shorter than these periods. This 
may result in that a collision in progress can be reversed at some conditions. Quantum na- 
ture of collisions is reflected in the formalism of quantum kinetics contrary to the Boltzmann 
kinetics dealing with free particle propagation between instant collisions. 



C. Physical meaning of time delays and advances 

In order to consider the variable shifts in the collision term in more detail, let us present 
the time delay as a sum of two quantities 

5t = ^Tdrift + 5TcoI = "Y - -fT^ (39) 

where Z = [vq — dKe T,^ / dp'^) ^ is the normalization factor which usually appears in de- 
scription of quasiparticle effects. In contrast to the quasiparticle case, here po and p are 
independent variables (not connected by any dispersion relation). The Bq function can be 
formally expressed as 

B„ = 2^. tani=-^. (40) 

The latter ratio is a measure of proximity of a virtual particle with given pq and p to the 
mass-shell, see Eq. (!30|) above. 

The value 5rdrift = -Bo/2 has the meaning of the drift delay time, 6tcoi = —Z^^/V, of the 
collision delay time, and 5r, of the effective total delay time. In the virial limit these time 



delays coincide with those studied in 
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201 ]. In particular, for quasiparticles Z > and 



M ^ 0, and therefore we obtain Bq = ZA = 2'kZ5{M) > and hence 5rdrift is positive 
delay, and St^oi is negative delay at arbitrary densities. To further clarify the meaning of 
these delays let us consider several examples. 



ttA^A system 
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If we consider propagation of a resonance, e.g., A isobar, in the virial limit of very 
small densities, 6 has the meaning of the irN phase shift. Thus in the virial limit, Bq of 
the resonance relates to the energy variation of the scattering phase shift (^3^) of the 
channel coupled to the resonance, cf. [igI. Due to the Galilean invariance the vacuum 
Green function and self-energy depend only on the energy in the ttN center-of-mass 
frame and Bq is recovered from the known dependence 6^^ (Ecm.) ■ As it follows from 
this dependence, i?o > and hence 5rdrift > (delay). There also appears a negative 
delay 6rcoi < (for Z > 0), since namely during the Z^^/T^ time the A resonance 
state forms the intermediate state of the vrA^ pair (causing an advance), compare with 
Ref. 2^. The 6t value in Eq. (|39l) has no definite sign. Its sign depends on which 
time is dominant, (5rdrift or 6tcoi- If we consider the pion propagation, Bq = Bq and in 
the virial limit 6 has the meaning of the in-medium A - nucleon hole phase shift. For 
the nucleon, Bq is related to ^ttA- 

Wigner resonances 

For a Wigner nonrelativistic resonance, when Re and F are independent of p, we 
have Z = 1. Then the value 6t in Eq. fl39l) becomes 



j£l-rv^ (41) 

which precisely coincides with the result Ref. |20!]. As it follows from f HTl) . in the 
vicinity of the resonance (M 0) 6t is positive (delay) and outside the resonance 
region (M^ ^ F^) 6t is negative (advance). If the interaction is attractive, the particle 
spends a shorter time in the interaction zone than that would be, if it moved with 
unchanged asymptotic velocity. This case corresponds to a delay. If the interaction 
is repulsive, the particle vise versa spends a longer time in the interaction zone, that 
corresponds to a negative delay. Thus, close to the Wigner resonance we may speak 
about an effective attraction and outside the resonance region, about an effective 
repulsion. 

Soft photon radiation 



Let us consider another example [2l|] of the photon radiation in the propagation of 
a source charge (say a proton) in neutral matter (e.g. the neutron one). The source 
particle looses the memory about collisions which have happened long ago (after a 
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time T ^ \Stcoi\ has past). The advance in the colhsion time can be understood in 
terms of multiple scatterings. The latter act coherently for r < |5tcoi| ~ 1/r. Indeed, 
for energies of the radiating quantum corresponding to the wavelength A-y = 2tt/uj. 



7 



> 



I^Tcoil the source particle simultaneously feels all rescatterings which it had and also will 
have on the time scale I^TcoiI. This results in the well known Landau- Pomeranchuk- 
Vligdal effect [22] experimentally proven at the Stanford Linear Accelerator Center 
23|. For large photon energies, is small and radiation occurs on a single center 
(neutron), whereas soft radiation occurs coherently on many centers. Thus appearance 
of the time advances 6tcoi in this example is associated with the multiple rescattering 
effects. 

In the above examples Z factors were assumed to be positive. However, Z can be negative 



24| . Change of the sign of 



in some energy-momentum regions far from the mass-shell, cf. 
Z is reflected in the respective sign change in 6tcoi- 

Moreover, interactions lead to the change in the generalized group velocity 

ttgroup = Z [v^ (42) 



dp 

that reflects in the time-space delays. To illustrate it let us consider an example of the pion 
propagation in the isospin-symmetric nuclear matter. For the sake of simplicity, we consider 
the pion propagation in a region close to the A-resonance and assume that the pion interacts 
with the nucleons only through the vrA^A-coupling /^ata- Then the real part of the pion 
self-energy becomes I2J] 

Here = tti^ — rriN, is the nucleon matter density, ao = const ~ 1. From here we can 
find Z^-factor. For uj —>■ uj^ we obtain 



Z-K 



2cja - 



J- A 



group 

^ 2kZ (44) 



resulting in an increase of the group velocity in the resonance region. Here it is worthwhile 
to note that for light in a dielectric medium there appears increase in the group velocity for 



packets with resonant frequencies [251] 
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D. Space and energy-momentum shifts 

As follows from Eq. (132|) the shift of the space variables is 

1 dRe 1 

5X = — — = 6Xdnit + SXcol = -B + Wgroup<5l"col- (45) 

We can further express the spatial shift fHSj) in terms of time delays (139|1 and velocity x of 
a test particle on its trajectory,'' cf. Eq. ([2] 



5x= X^Tdrift + 'Ugroup^T-col. (46) 

Eq. fj46l) demonstrates that to a delayed collision the test particle moves along its trajectory. 
Therefore, the drift time delay 5rdrift unambiguously results in a definite space shift x^Tdrift- 
The collision itself is associated with an additional time delay ^Tcoi, which implies that 
the collision is not instant, as it is treated in the kinetic equation, but requires certain 
time for complete decoupling from intermediate states (e.g., the pion spends some time in 
the intermediate A-nucleon-hole state, a soft photon requires certain time to be formed in 
multiple collisions of the proton with neutrons). Therefore, this additional delay gives rise 
to an additional shift 'Ugroup'^Tcoi of the particle with respect to its "collisionless" trajectory 

(EH). 

Since during time-delay interval 5t a particle propagates in the mean field Re S^, it also 
changes its energy and momentum. The energy-momentum shift fl32l) can be represented in 
a similar form 

5p^^ = p'^^Tdrift - Z—^5t,,,, (47) 

where = {E, p) is the rate of change of the 4-momentum of a test particle on its trajectory, 
cf. Eqs. ( 128|) and (!29|) . We again distinctly see two contributions to The first one is 
associated with the energy-momentum change during the path of the particle to a delayed 
collision. The second one is related to an intrinsic nonlocal character of the collision itself. 
Thus, we see that all shifts in space and energy-momentum variables are consistent with 
time delays discussed above. 

The energy shift can be expressed in the alternative form 

M2 -rV4 fdtReJ:^ 



^ Here for briefness we omit the unnecessary subscript i of the test particle. 
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where we used Eq. (|9]). In the case of the Wigner resonance, the first term on the r.h.s. of 
this relation is —5t ■ c^T-Re S^. Thus, the energy-shift scale is ~ |Re S-^ ■ 5r|/rch, where Teh is 
the characteristic time of the self-energy variation. In the virial limit ReS^ = 47rpReF(0), 
where -F(O) is the forward scattering amplitude in the energetic normalization, p is the 
density of the medium. Then we obtain 

6pQ 47rRe F(0)5r ■ dtp, (49) 

thus relating the value Teh to the variation of the density. In the kinetic region the scale 
is the collisional time, r^h ~ |5tcoi|, and hence the correction 5pQ can be not small. In the 
region of the validity of the hydrodynamics the density varies on a much longer scale than 
the collisional one, and hence the correction 5pQ is small but still does not vanish. When 
the global equilibrium is reached this correction disappears. 



Since the time delay occurs in all the quantities entering the collision term, and the 
value of the 5t delay is finite even in the case of global equilibrium (e.g., see Eq. (HT!) ). 
the thermodynamic quantities, such as the pressure, energy, etc., also acquire corrections 
proportional to this time delay. 

At the general level, the relation of thermodynamic properties to the time delays can be 
understood in terms of ergodicity. The system spends certain time in a certain phase-space 
region, which is proportional to the density of states in that region. Change in this time cor- 
responds to the change in the density of states and thereby to the change in thermodynamic 
properties, cf. |20|. We have generalized this statement to the unified 8-phase-space. 

Another important feature of the kinetic description is approaching to thermal equilib- 
rium during the evolution of a closed system. In terms of transport theory it means that the 
H-theorem takes place for transport equations, which implies construction of the entropy 
related to the kinetic equation. These problems have been considered both for the BM 
equation and for the KB equation , [lo| . Since the NL form of the quantum kinetic 

equation (|3T|) is equivalent to the KB equation including second-order terms in space-time 
gradients, the NL entropy differs from the KB entropy only in the second-order gradient 
terms. This reduces the problem of the H-theorem for the new NL kinetic equation to the 

n 

already solved problem for the KB equation |10l |. 
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IV. SUMMARY 

We have considered properties (advantages and disadvantages) of available formulations 
of the off-shell kinetic equation, i.e. those in the Kadanoff-Baym and Botermans-Malfliet 
forms, as well as a new nonlocal form which we derived in this paper. Within the range of 
formal applicability, all these forms are equivalent to each other. Under this applicability 
range we mean states of a system which are close to equilibrium, where the space-time 
gradient expansion is strictly applicable. In particular, it implies that conservation laws of 
Noether currents and energy-momentum are satisfied at least to the extent of validity of the 
dynamics, i.e. up to zero-order gradient^. Note that these conservation laws are even exactly 
fulfilled for the Kadanoff-Baym form [9]. However, these equations are frequently applied 
beyond the scope of their applicability, e.g., when the collision term is large. Such a situation 
takes place, e.g., in description of the initial stage of heavy-ion collisions. In such situations 
we certainly wish to respect conservation laws which allow us to keep control of numerical 
codes. From this point of view the KB form is certainly preferable. As the conservations 
are exact, we can still use the gradient approximation, relying on a minor role of this rather 
short initial stage of heavy ion collisions in the total evolution of a system. However, the 
efficient test-particle method is not applicable for solution of the Kadanoff-Baym equation. 

The nonlocal kinetic equation differs from the Kadanoff-Baym one only in third-order gra- 
dient terms, while the Botermans-Malfliet form, in the second order. Due to this the nonlocal 
form of the kinetic equation is certainly preferable as compared with the Botermans-Malfliet 
one, since it more accurately conserves the Noether currents and the energy-momentum. 
Deviations from the Noether quantities appear only in the second gradient order, as com- 
pared with flrst order for the Botermans-Malfliet equation. In addition, the nonlocal form 
allows application of the test-particle method very similarly to that has been done for the 



Botermans-Malfliet form 
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14| . Presenth 



;here are numeric schemes dealing with delayed 



(advanced) collisions of the test particles [20 



261 ]. Thus, the nonlocal form of the quantum 



kinetic equation is a reasonable compromise between the Kadanoff-Baym and Botermans- 
Malfliet equations, which is practical like the Botermans-Malfliet form and at the same time 
reasonably preserves conservations of Noether quantities. 

Physical meaning of delays and advances in 8-dimensional phase space appearing in the 
nonlocal form of the kinetic equation is clarifled. It is demonstrated that in addition to 
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the practical utility the nonlocal form provides us with clear physical interpretation of var- 
ious terms in the Kadanoff-Baym equation and also allows us to consider delicate physical 
phenomena, which are beyond the scope of the Botermans-Malfliet approximation. 
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Appendix A: MATRIX NOTATION 

In calculations that apply the Wigner transformations, it is necessary to decompose 
the full contour into its two branches — the time-ordered and anti-time- ordered branches. 
One then has to distinguish between the physical space-time coordinates x, . . . and the 
corresponding contour coordinates which for a given x take two values x" = {x~) and 
x'^ — {x^) {fj. G {0, 1, 2, 3}) on the two branches of the contour (see figure 1). 

t: 



y 

-•- 

t 



to t+ oo 

Figure 1: Closed real-time contour with two external points x,y on the contour. 

Closed real-time contour integrations can then be decomposed as 

r poo pto poo poo 

dx . . . = dx . . . -\- dx . . . — dx . . . — dx . . . , 

J C J to J oo J to J to 



(Al) 



where only the time limits are explicitly given. The extra minus sign of the anti-time-ordered 
branch can conveniently be formulated by a { — h} "metric" with the metric tensor in { — h} 
indices 

= K) = (J_°J (A2) 
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which provides a proper matrix algebra for multi-point functions on the contour with "co"- 
and "contra" -contour values. Thus, for any two-point function F, the contour values are 
defined as 



Fij{x,y) :-- 



F{x\y^), i,je {-,+}, with 
a,kF''^{x,y), F\.(x,y) := F'\x,y)auj 



(Tik(J.jiF^\x,y), (j'y = 6, 



ik 



(A3) 



on the different branches of the contour. Here summation over repeated indices is implied. 
Then contour folding of contour two-point functions, e.g., in Dyson equations, simply be- 
comes 



H{x\y'') = H'\x,y) = [F G] 
= j dzF{x\z)G{z,y^ 



ik 



dzF\{x,z)G^\z,y) 



(A4) 



in the matrix notation. 

Due to the change of operator ordering, genuine multi-point functions are, in general, 
discontinuous, when ever two contour coordinates become identical. In particular, two-point 
functions like iF{x,y) = (TcA{x) B {y) \ become^ 



iF{x,y) 



^iF-{x,y) iF-+{x,y)^ 
yiF+-(x,?/) iF++(a;,?/) J 



( (TA{x)B{y)) T{B{y)A{x] 



\ 



A{x)B{y) 



T-'A{x)B{y) 



\ 



J 



(A5) 



where T and are the usual time and anti-time ordering operators. Since there are 
altogether only two possible orderings of the two operators, in fact given by the Wightman 
functions F~^ and F^~ , which are both continuous, not all four components of F are 
independent. Eq. (lASp implies the following relations between nonequilibrium and usual 
retarded and advanced functions 



F^{x,y) 
F^{x,y) 



F-ix, y) - y) = y) - y) 

0(xo - yo) {F^'ix.y) - F-+{x,y)) , 

y) - y) = F~+{x, y) - y) 

-Q{yo-xo) {F+-{x,y) - F~+{x,y)) , 



(A6) 



Frequently used alternative notation is F< — F ^ and — 
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where Q{xq — yo) is the step function of the time difference. The rules for the co-contour 

functions F etc. follow from Eq. flA3l) . 

For slightly inhomogeneous and slowly evolving systems, the degrees of freedom can be 
subdivided into rapid and slow ones. Any kinetic approximation is essentially based on 
this assumption. Then for any two-point function F'^^{x,y), one separates the variable 
^ = (ti — t2,ri — r^), which relates to rapid and short-ranged microscopic processes, and 
the variable X = ^{ti +t2, + ^2), which refers to slow and long-ranged collective motions. 
The Wigner transformation, i.e. the Fourier transformation in four-space difference ^ = x—y 
to four-momentum p 

r^{X;p) = I d^e'P^r^X + ^/2,X-a2), t,j G {-+} (AT) 

where X = (t, x), leads to the corresponding Wigner densities in four-phase-space. The 
gradient expansion converts the Wigner transformation of any convolution of two-point 
functions into a product of the corresponding Wigner functions plus higher order gradient 
terms 

Jd^e'P^ydzfix,zMz,y)^ = (^exp '-^ {d,dx' - dxd,,) f{X,pMX',p')^ ^ JA8) 

^ fiX,pMX,p) + '-^{f{X,p),^{X,p)}, (A9) 
where the first order terms are given by Poisson brackets 

,M0) 

here in covariant notation. We would like to stress that the smallness of the hdx ■ dp comes 
solely from the smallness of space-time gradients dx, while momentum derivatives dp are 
not assumed to be small. This point is sometimes incorrectly treated in the literature. 
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